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ABSTRACT 

An extratropical cyclone that crossed the United States on 9-11 April 2009 was successfully simulated at 
high resolution (3-km horizontal grid spacing) using the Colorado State University Regional Atmospheric 
Modeling System. The sensitivity of the associated warm front to increasing pollution levels was then explored 
by conducting the same experiment with three different background profiles of cloud-nucleating aerosol 
concentration. To the authors" knowledge, no study has examined the indirect effects of aerosols on warm 
fronts. The budgets of ice, cloud water, and rain in the simulation with the lowest aerosol concentrations were 
examined. The ice mass was found to be produced in equal amounts through vapor deposition and riming, and 
the melting of ice produced approximately 75% of the total rain. Conversion of cloud water to rain accounted 
for the other 25%. When cloud-nucleating aerosol concentrations were increased, significant changes were 
seen in the budget terms, but total precipitation remained relatively constant. Vapor deposition onto ice 
increased, but riming of cloud water decreased such that there was only a small change in the total ice pro- 
duction and hence there was no significant change in melting. These responses can be understood in terms of 
a buffering effect in which smaller cloud droplets in the mixed-phase region lead to both an enhanced vapor 
deposition and decreased riming efficiency with increasing aerosol concentrations. Overall, while large 
changes were seen in the microphysical structure of the frontal cloud, cloud-nucleating aerosols had little 
impact on the precipitation production of the warm front. 


1. Introduction 

Extratropical cyclones are the primary transporters of 
moisture, heat, and energy poleward in the midlatitudes 
(e.g., Trenberth and Stepaniak 2003; Schneider et al. 
2006). The fronts associated with extratropical cyclones 
contribute approximately 68% of the total precipitation 
at midlatitudes (Catto et al. 2012). Additionally, their 
deep and extensive cloud shields have large impacts for 
radiation locally and on climate scales (e.g., Ramanathan 
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et al. 1989; Stewart et al. 1998). Should aerosols be able 
to alter the radiative, microphysical, and/or dynamical 
properties of extratropical cyclones, they could have 
important consequences for global change scenarios. 
However, very few studies have investigated the indirect 
effect of aerosols on extratropical cyclones, despite the 
fact that both changes in pollution emissions globally 
and local seasonal enhancements of aerosol concentra- 
tions because of burning could impact these systems, as 
shown conceptually by Wang et al. (2009). 

Aerosols can affect cloud systems by changing the 
microphysical structure of clouds. A change in cloud 
microphysical structure can affect the radiative proper- 
ties of clouds (Twomey 1977), change the lifetime of 
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clouds (Albrecht 1989), or alter the rates of microphysical 
processes associated with precipitation production and 
latent heat release. Changes in diabatic heating rates can 
then induce changes in vertical velocities and other dy- 
namical fields. Understanding all of these indirect effects 
and feedbacks of aerosols is an ongoing and difficult task, 
especially since the environment and dynamics of the 
system itself can alter the response to aerosols (Matsui 
et al. 2006; Fan et al. 2009; Khain 2009; Storer et al. 
2010). Since the mesoscale dynamics of an extratropical 
cyclone differ greatly between the cold and warm fronts, 
a detailed study of all of the impacts of aerosols on 
extratropical cyclones is too extensive for one study; 
therefore, this paper will focus on just one part of the 
extratropical cyclone, the warm front. 

While there are no studies reported in the literature 
that specifically examine the influence of aerosols on 
warm fronts, aerosol impacts on other mixed-phase 
cloud systems have been investigated, and the results of 
these studies are likely to be relevant to this study. 
Cheng et al. (2010) examined a mixed-phase front sys- 
tem in Taiwan and found that precipitation could be 
either increased or decreased with increasing cloud 
condensation nuclei (CCN) depending on changes in the 
depositional growth of ice and riming. CCN in mixed- 
phase clouds have been shown to result in decreased 
riming efficiency in observations (Borys et al. 2003) and 
modeling studies of mixed-phase orographic clouds 
(Saleeby et al. 2009; Xue et al. 2012) because of de- 
creased supercooled cloud droplets sizes. Both obser- 
vations and modeling studies of these clouds show 
changes in the spatial distribution of precipitation (e.g., 
Givati and Rosenfeld 2004; Jirak and Cotton 2006; Lynn 
et al. 2007; Saleeby et al. 2009) but do not always agree 
on the magnitude or sign of the total precipitation 
change, which can depend on environmental factors and 
the extent of the mixed-phase region (Muhlbauer et al. 
2010; Zubler et al. 2011). Many results of mixed-phase 
convective cloud studies indicate that higher CCN con- 
centrations can cause more liquid water to be lofted 
above the freezing level leading to additional freez- 
ing and latent heat release (Khain et al. 2005; van den 
Heever et al. 2006). The enhanced buoyancy and feed- 
backs to the updraft are often referred to as the aerosol 
invigoration effect, but this effect has not been studied 
extensively in nonconvective cloud types. 

In addition, should aerosols be able to significantly 
alter latent heating rates, results of studies that examine 
the sensitivity of warm fronts to latent heating may also 
be relevant to this study. First, modeling studies have 
shown that latent heat processes, including condensa- 
tion and melting, lead to mesoscale banded features 
over the warm-frontal surface (Hsie et al. 1984; Szeto 


and Stewart 1997). Evaporational and sublimational 
cooling can enhance regions of mesoscale descent (Huang 
and Emanuel 1991; Clough et al. 2000). Through these 
modifications of vertical velocity and circulation, as well 
as through local changes in the thermal structure, latent 
heating has been seen to accelerate frontogenesis and 
frontal propagation speed for both warm and cold fronts 
(Hsie et al. 1984; Szeto and Stewart 1997; Reeves and 
Lackmann 2004). 

Before the impacts of aerosols on warm fronts can be 
fully understood, the budgets of ice and liquid hydro- 
meteor mass in warm-frontal clouds must be examined 
since aerosol indirect effects vary with cloud phase. Such 
budgets of cloud ice and rain, but not cloud water (here 
and throughout the study, “cloud water” means “cloud 
liquid water”), have been studied before by Rutledge 
and Hobbs (1983) and Gedzelman and Arnold (1993) 
for warm fronts; however, both of those studies used 
idealized, two-dimensional model setups with more 
primitive microphysical schemes than is used in the 
present study. The current work will expand on these 
previous studies by presenting budgets of ice, rain, and 
for the first time cloud water for warm-frontal clouds 
from a three-dimensional, high-resolution cloud-resolving 
simulation that employs sophisticated microphysics 
and aerosol parameterization schemes. The sensitivity 
of the budgets and latent heating profiles of the warm 
front will then be tested by altering the background 
profile of cloud-nucleating aerosol concentration in two 
additional simulations. 

In section 2, the model will be described and a syn- 
optic overview of the selected storm given. In section 3, 
data analysis methods will be described. Section 4 in- 
cludes the budgets of ice, cloud water, and rain, in warm- 
frontal clouds. In section 5 the sensitivity of these budgets 
and the warm front as a whole to perturbations in aerosol 
concentrations is discussed. 


2. Case overview 

a. Case selection and model setup 

To approach this problem, an extratropical cyclone 
was simulated using the Regional Atmospheric Modeling 
System (RAMS) (Pielke et al. 1992; Cotton et al. 2003). 
The case selected was the extratropical cyclone that de- 
veloped in the lee of the Rocky Mountains and tracked 
across the United States from 9 to 11 April 2009. This case 
was selected for several reasons. Most importantly, it had 
a long-lived and well-defined warm front. Additionally, 
no other extratropical cyclone was near enough to in- 
teract with this storm, which would further complicate the 
analysis. Third, the case was chosen because it did not 
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move from water to land or vice versa until after the 
simulation period. In this way transitions between water 
and land surfaces do not complicate interpretation of the 
results. 

The simulations were set up with an outer and inner 
two-way nested grid, the locations of which are shown in 
Fig. la. The outer grid has 250 X 300 grid points with 
15-km horizontal spacing (3750 km X 4500 km) while the 
inner grid has 697 X 1027 grid points with 3-km hori- 
zontal spacing (2091 km X 3081 km). While portions of 
fronts and extratropical cyclones have been run at finer 
resolution (e.g., Garvert et al. 2005; Molthan and Colie 
2012), to our knowledge, no extratropical cyclone in its 
entirety has been simulated at such a high grid resolution 
before. Both grids have 45 vertical sigma levels, with 75-m 
spacing near the surface stretching to 1-km spacing at the 
model top. The time steps on the inner and outer grids 
were 5 and 20 s, respectively. A cumulus parameteriza- 
tion scheme based on Kuo (1974) and Molinari (1985) 
was implemented on the outer grid only. The Harrington 
(1997) radiation scheme and a Smagorinsky (1963) tur- 
bulence scheme with modifications by Lilly (1962) and 
Hill (1974) were implemented on both grids. 

The simulations use a two-moment bin-emulating mi- 
crophysical scheme that includes five prognostic ice spe- 
cies: pristine ice, snow, aggregates, graupel. and hail, and 
three prognostic liquid species: two cloud droplet modes 
and rain (Meyers et al. 1997; Saleeby and Cotton 2004, 
2008). Graupel and hail are distinguished by density rather 
than size in the RAMS microphysics scheme (Walko et al. 
1995). Their densities are defined to be 0.3 and 0.9 g cm~ 3 , 
respectively. These densities for graupel and hail allow 
for a transition in density as ice particles become more 
heavily rimed or gradually melt and therefore allow 
the fall speeds of these particles to be more accurately 
computed. Graupel and hail both have diameters of 
approximately 1-2 mm in these simulations; however, 
the hail concentration is about 4 times larger than that of 
graupel and dominates the mixing ratio because of its 
higher density. While hail is not typically associated with 
warm fronts, the “hail” produced by the model is more 
appropriately described as heavily rimed ice or partially 
melted ice, particularly below the melting level. Both 
graupel and hail can be, but are not necessarily, partially 
melted. Therefore, we will refer to them as the “mixed- 
phase” species throughout the remainder of the text. 

By using a bin-emulating scheme, each species is 
binned by size for the calculation of collection effi- 
ciencies and fall speeds; these calculations give a more 
accurate representation of riming and autoconversion 
than can be simulated by using a traditional bulk mi- 
crophysics scheme. This parameterization scheme is 
therefore well suited for the study of the microphysics 


and precipitation of the warm front without being too 
memory intensive. 

RAMS explicitly tracks the rates of many cloud pro- 
cesses. Those that will be presented here (the names 
used in the figures are given in parentheses) are total 
vapor deposition and sublimation of ice species (vapor 
to ice), total condensation and evaporation of cloud 
species (vapor to cloud) and rain (vapor to rain), cloud 
droplet nucleation (cloud nucleation), conversion of 
cloud water to rain through autoconversion of cloud 
droplets and accretion of cloud droplets by rain (cloud to 
rain), melting (melting), riming of cloud droplets (rim- 
ing of cloud) and raindrops (riming of rain), and col- 
lection of ice by rain resulting in the melting of the ice 
(ice to rain). These allow for the budgets of ice, cloud 
water, and rain to be explicitly calculated in the simu- 
lations. Note though that the two processes comprising 
each of the following pairs — deposition and sublimation, 
condensation and evaporation, and the cloud to rain 
processes — are not tracked separately and thus cannot 
be examined individually. 

The model was initialized using Global Forecast System 
(GFS) analysis at 0000 UTC 9 April 2009 and integrated 
for 48 h. The lateral boundaries were nudged using ad- 
ditional GFS analyses. The profile of those aerosols that 
can serve as cloud condensation nuclei (CCN) was maxi- 
mized at the surface and decreased linearly from the 
surface value up to 4 km above which it was constant at 
100 cm~ 3 . The model was run 3 times with surface aerosol 
concentrations of 400, 800, and 1600 cm~ 3 ; these con- 
centrations are representative of clean, average, and pol- 
luted conditions over continents (Andreae 2009). These 
sensitivity runs will be referred to as CCN400, CCN800, 
and CCN1600, respectively. Apart from the differences in 
initial aerosol concentrations, the model setups were 
identical. The cloud-nucleating aerosol particles were al- 
lowed to be advected horizontally and vertically but were 
not depleted during cloud nucleation. To prevent runaway 
nucleation, aerosol particles were nucleated under ap- 
propriate environmental conditions only if their number 
concentration did not exceed the number concentration of 
cloud droplets present in the same grid box. Had the 
cloud-nucleating aerosol been depleted, most particles 
would have been rained out quickly, leaving a very clean 
environment in all three cases since no surface sources of 
aerosol were present. Additionally the aerosols were not 
radiatively active, so no assessment of the aerosol direct 
effect is made in this study. 

b. Synoptic overview 

Figure 1 shows output from CCN400 depicting the 
evolution of the storm. At 1200 UTC 9 April, a broad 
trough exists at both 300 and 500 hPa over the Rocky 
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Fig. 1. (a),(c),(e),(g) Plots of 300-hPa wind speed (shaded; ms 1 ) and 500-hPa geopotential height (contoured; m) 
and (b),(d),(f),(h) surface potential temperature (shaded; K) and geopotential height (contoured; m) output from the 
outer grid at (a),(b) 1200 UTC 9 Apr, (c),(d) 0000 UTC 10 Apr, (e),(f) 1200 UTC 10 Apr, and (g),(h) 0000 UTC 
11 Apr. The red and brown lines in (b),(d) are warm fronts and drylines, respectively, and the blue line in (f) is a cold 
front analyzed by hand. The thick pink line in (f) is the objective location of the warm front. Fronts emanating from 
the storm center were not detected at the final time shown in (h). 
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Mountains associated with a closed low at 850 hPa be- 
ginning to form in eastern Colorado and western Kan- 
sas. The trough at 500 hPa has formed a closed low, the 
downstream ridge has weakened, and the low at 850 hPa 
has deepened 12 h later. There is an enhanced potential 
temperature gradient at 850 hPa extending eastward 
from the low center, consistent with the presence of 
confluent background flow (Schultz et al. 1998). This is 
the warm front that will be the focus of this study. There 
is a tongue of warm air reaching the center of the cir- 
culation from the south and west that is associated with 
a dryline (Figs. lb,d). This boundary was associated with 
a severe weather outbreak in the southeast on 9-10 
April but will not be discussed in this paper. 

At 1200 UTC 10 April the system is near its peak in- 
tensity. The low centers have become vertically stacked 
and begin to weaken shortly after this time. The warm 
front of interest is now associated with a distinct trough 
in the 850-hPa field and its eastern end is beginning to 
interact with the Appalachian Mountains. The jet streak 
at 300 hPa in the northeast has intensified but shifted 
little in the horizontal and is likely aiding in large-scale 
lift over the warm front. By 0000 UTC 11 April the low 
center is just to the west of the mountains. As it is not the 
objective of this study to examine the influence of the 
mountains on the warm front, the simulation was ended 
here, though the system did continue to the Atlantic 
Ocean and regained strength over the Gulf Stream. 

Figure 2 shows a time- and space-averaged, front-normal 
cross section of cloud ice, cloud water, and rain from the 
CCN400 run of the warm front. Details on how these 
cross sections were created are described in section 3. 
The isentropes indicate that the warm front is shallow 
with a slope of about 1:150 and that the frontal cloud is 
predominantly stratiform. It can be seen that the mixed- 
phase region is very deep, extending from about 2 to 
6 km and that the melting level located near 2-3 km 
throughout most of the system. 

c. Validation 

Although we are simulating an actual extratropical 
cyclone, the goal of this research is not to present an in- 
depth analysis of a case study. Rather, we are interested 
in the sensitivity of the warm front to perturbations in 
cloud-nucleating aerosol concentrations. As a conse- 
quence, the results are treated more like those of an 
idealized modeling study, where the case study data are 
simply used to simulate the storm type of interest, in this 
case an extratropical cyclone. Therefore, it is not im- 
portant that the model accurately capture the time 
evolution or exact structure of the cyclone, only that the 
model simulates an extratropical cyclone that shows 
average characteristics that are representative of such 


Cloud Ice 



Distance from Front (km) 



Fig. 2. Time-averaged cross-sections of (a) cloud ice, (b) cloud 
water, and (c) rain mixing ratio (g kg -1 ). Thin black lines are 
isentropes (K), and the thick black line is the 0°C line or melting 
level. The gray, pink, and red lines show the 0.005 g kg -1 contours 
of cloud ice, cloud water, and rain, respectively, for reference. 
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FlG. 3. (a) Radar-derived estimate of the 24-h precipitation total ending 1200 UTC 10 Apr from the National 
Weather Service, (b) Precipitation total for the same time period, but from RAMS output. 


a system. However, it is useful to know that the model 
physics are to a first order correct and that the model has 
captured some basic features of the actual system. Total 
precipitation is a good field to use for this since it con- 
tains information about the storm location as well as 
rainfall. 

Figure 3a shows the National Weather Service (NWS) 
Multisensor Precipitation Estimator (Lawrence et al. 
2003) rainfall and Fig. 3b the CCN400 total precipitation 
for the 24-h period ending 1200 UTC 10 April on the 
inner grid. The general pattern of precipitation is well 
captured, especially around the warm front. RAMS cor- 
rectly places the maximum in precipitation in northern 
Missouri and southern Illinois and Indiana, though per- 
haps it is shifted slightly north and east. Comparisons of 
model output to NWS surface analyses also show a shift 
to the northeast of the modeled location of the warm 
front and cyclone center relative to their observed loca- 
tions (not shown). The local maxima in precipitation in 
RAMS are higher than those seen in the radar estimate, 
but the average totals of 1-4 cm seem to be consistent. 
RAMS does underpredict rainfall in the southeast asso- 
ciated with the secondary warm front, but that is not 
a major concern for the analysis that will be presented 
here as our focus is on the primary warm front. 

3. Methods 

a. Warm front detection 

To objectively locate the warm front at each hour, we 
used modeled potential temperature 0 surf and geostrophic 
wind velocity V r , surf, and applied the front detection al- 
gorithm of Hewson (1998) on equal-sigma surfaces rather 
than equal-pressure surfaces as in previous work (e.g., 
Naud et al. 2012). While Hewson (1998) recommended 


using the 900-hPa level for detecting surface fronts, we 
found that the surface fields were much less noisy and 
produced cleaner results than other near-surface levels. 
The technique proposed by Hewson (1998) was originally 
conceived for fields of about 100-km horizontal resolu- 
tion, but here the resolution is 3 km, so we had to make 
several adjustments. We first subsampled 0 surf and UG_surf 
every 30 km and then applied a smoothing routine. These 
two fields were then fed into the front detection routine. 
Only those fronts identified by the routine with potential 
temperature gradients greater than 25 X lO^KnU 1 
were kept, and the sign of the product of geostrophic wind 
velocity and potential temperature gradient delineated 
between warm and cold fronts. 

Despite the upscaling and smoothing, multiple warm 
fronts were detected, and additional masks were applied to 
objectively isolate the primary front. Warm-frontal points 
along coastlines were removed. These points existed sim- 
ply because of the large contrast in potential temperature 
between land and ocean and did not represent the warm 
front associated with the extratropical cyclone. Warm- 
frontal points were also required to have cloud cover 
within 60 km to the north in order to eliminate the 
boundary oriented north-south in the warm sector of 
the storm (Figs. lb,d) that was associated with the se- 
vere weather outbreak mentioned earlier. Next, warm- 
frontal points with four or fewer neighbors within 
150 km, one or zero neighbors within 60 km, or zero 
neighbors within 300 km in the next hour were thrown 
out. The effects of these conditions were to remove 
short-lived warm-frontal structures. Finally, a 10° poly- 
nomial was fit to the remaining warm-frontal points in 
order to preserve the shape and location of the warm 
front while smoothing clusters of points and to recover 
the original 3-km inner grid resolution. 
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Fig. 4. Hovmoller diagrams of the model average u, v, and potential temperature at the surface as a function of distance from the warm 

front. See text for further details. 


The final objective locations of the warm front at each 
hour (see, for example, the one in Fig. If) show that the 
detection algorithm effectively locates the warm front 
such that it is continuous in time. It performs best from 
hours 27-42 of the simulation during which time almost 
the entire length of the front propagates forward at each 
time step. It is over this interval that the analyses pre- 
sented here will be conducted. 

b. Averaging 

The objective warm-front locations were used to create 
average cross sections through the front. The distance to 
the storm center and length of the front varies in time. 
However, to maintain a consistent data volume for each 
time, warm fronts were defined to begin 60 km east of the 
low pressure center (the largest distance between the low 
pressure center and the first point of any objective front) 
and to be 573 km long (the length of the shortest front in 
the simulation). The warm fronts are nearly parallel to 
the west-east direction at each time step (Fig. If, for ex- 
ample), so no rotation of the data was necessary. Model 
output was taken along a north-south transect at each 
grid point along the front and averaged along the east- 
west direction to obtain a single composite cross section 
perpendicular to the front at each hour. The meridional 
span captured the bulk of the precipitation and cloud field 
associated with the warm front. Figure 4 shows Hovmoller 
diagrams of the zonal wind u, the meridional wind v, and 
potential temperature at the surface as a function of 
distance from the front. It can be seen that the warm- 
front detection and averaging correctly places the front 
on the warm side of the baroclinic zone and that the 
cross-front wind shifts from positive to negative as ex- 
pected. This indicates that the front detection and av- 
eraging is effectively capturing the features of the front. 

In many studies of aerosol indirect effects, quantities 
are averaged only over regions where some definition 
of cloud is met. However, in this study no minimum 
threshold for condensate is used in any of the averaging 
of cloud processes or cloud properties since cloud and/or 


precipitation are present in most of the subsetted do- 
main (Fig. 2). In this way, the number of points being 
averaged is the same in each case and any differences 
between simulations must be due to actual changes in 
the quantity being shown, rather than changes in sample 
size. The exceptions are hydrometeor diameters, which 
are averaged only where they are nonzero and are 
weighted by the number concentration of the proper 
species. All vertical profiles shown in the figures to fol- 
low are averaged from the surface location of the front 
to 525 km north of the front. Areas immediately south 
of the surface front did not receive substantial pre- 
cipitation and therefore were not included in the aver- 
ages. By averaging to 525 km ahead of the front, nearly 
the full extent of the warm-frontal cloud is captured, 
though at its farthest reaches, only ice exists and liquid 
species are not present (Fig. 2). 

4. Cloud budgets 

Figure 5 shows average vertical profiles of each of the 
most important processes contributing to the production 
or depletion of ice species, cloud water, and rain in 
CCN400. The ice budget (Fig. 5a) shows that vapor 
deposition is the dominant growth process above about 
4.5 km whereas riming is dominant below that level. 
This being a warm spring case (Fig. 1), very little ice 
reaches the surface in the simulation and therefore most 
ice that is produced must eventually melt. Sublimation 
and the collection of ice by rain are minor sinks of ice 
compared to melting. These results lie between the two 
cases examined by Rutledge and Hobbs (1983). In the 
feeder zone of their weak updraft case they found no 
contribution to the ice budget by riming, whereas in the 
feeder zone of their mesoscale updraft case they found 
the riming rate to be 2-3 times greater than the de- 
positional growth rate. 

Neither the study by Rutledge and Hobbs (1983) nor 
that by Gedzelman and Arnold (1993) examined the 
cloud water budget. The cloud water budget in this study 
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Fig. 5. Vertical profiles of each of the major processes contributing to the budget of (a) ice, (b) cloud water, and (c) rain. Note the melting 
level is near 2 km. Positive (negative) values indicate growth (decay) of the species. See the text for a full description of each process. 


(Fig. 5b) shows two peaks in production, one in the 
mixed-phase region of the cloud near 4 km and one just 
below the melting level near 2 km. The upper peak ap- 
pears to be caused by updrafts creating supersaturated 
conditions while the lower peak may exist because of 
supersaturated conditions being created by the latent 
cooling due to melting. The cloud nucleation rate is 
slightly greater than the condensational growth rate at 
both peaks and significantly greater in the region be- 
tween the peaks. The main process contributing to the 
depletion of cloud water at each peak is different. At the 
upper peak, cloud-droplet-to-rain conversion rates are 
small and riming rates dominate. However at the lower 
peak, cloud droplet to rain conversion rates are larger 
than the riming rate, though riming is still significant at 
the melting layer. Cloud droplet sizes and rain mixing 
ratios are both much smaller at the upper peak (Fig. 6), 
which would reduce the rates of autoconversion and 
accretion. When integrated through the column, riming 
contributes about 33% more to the depiction of cloud 
water than do cloud to rain conversion processes. 

Finally the rain budget is shown in Fig. 5c. Melting of 
ice is by far the most important producer of rainwater, 
a result found by both Rutledge and Hobbs (1983) and 
Gedzelman and Arnold (1993). The conversion of cloud 
water to rain is of secondary importance. The budgets 
show that about 75% of rain mass is due to melting of ice, 
where half of the ice mass is formed through depositional 
growth and half through the riming of cloud water. The 
remaining roughly 25% of rain mass is produced through 
warm rain processes in this case. The seeder-feeder 
process is thought to be responsible for the enhancement 


of precipitation in warm-frontal rainbands (Houze et al. 
1981; Rutledge and Hobbs 1983). Though no distinction 
has been made between banded and nonbanded pre- 
cipitation in this study, these results also suggest that ice 
crystals falling through the mixed-phase region of the 
cloud efficiently collect cloud water that would not oth- 
erwise be converted to rain; this collection contributes 
significantly to the total precipitation. 

5. Sensitivity to changes in cloud-nucleating 
aerosol concentrations 

a. Microphysical cloud properties 

The sensitivity of the mass budgets and cloud proper- 
ties to increases in cloud-nucleating aerosol concentra- 
tions will now be discussed by comparing results from the 
CCN400, CCN800, and CCN1600 simulations. Figure 6 
shows average vertical profiles of hydrometeor number 
concentration, mixing ratio, and diameter in the top three 
rows. The percentage change from CCN400 to CCN1600 
in these quantities is shown in the bottom row. As dis- 
cussed in section 2a, pristine ice, snow, and aggregates 
(PSA) have been averaged together, the mixed-phase 
species (GH) are combined, as are the two cloud droplet 
modes. Although the responses of each ice species to the 
sensitivity tests differ, grouping them in this way shows 
the overall impact of changes in aerosol concentrations. 
The PSA number concentration, diameter, and mixing 
ratio are dominated by aggregates below 8 km. It is 
important to note that profiles of mixing ratio are not 
expected to match the total budget profiles in Fig. 5 
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Fig. 6. Average vertical profiles of (top row) number concentration, (second row) mixing ratio, and (third row) number concentration 
for each of the four groups of hydrometeors described in the text. Diameters were averaged only where they were nonzero and weighted 
by number concentration, but no minimum threshold was used in averaging number concentrations or mixing ratios, (bottom row) The 
percentage change of CCN1600 from CCN400 for number concentration, mixing ratio, and diameter of each hydrometeor group. 


since advection and sedimentation were not included in 
the budget calculations. 

With increasing pollution cloud droplets are more nu- 
merous and smaller which is expected according to the 
cloud albedo effect (Twomey 1977), and the overall cloud 
water path is increased. PSA mixing ratio is highest in 
CCN1600 below 8 km, though CCN400 and CCN800 are 
similar; the trend becomes clearly monotonic at 4 km and 
below. Fewer GFI particles exist with enhanced aerosol 
concentrations although the differences in number are 


nearly zero below the melting level. Overall there is less 
mass in this category above the melting level in the more 
polluted cases, but again, below the melting level dif- 
ferences are very small. Throughout the column, rain 
exhibits very little variation in mixing ratio among the 
three simulations. There are monotonic trends in number 
concentration above 2 km where drops are less numerous 
but below 2 km, the trend is nonmonotonic. There is also 
no appreciable variation in rain droplet size in the lowest 
4 km of the column. 
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Fig. 7. Vertically integrated and averaged rates of each of the major water budget processes. 
The two numbers printed with each process are the percentage change of CCN800 and 
CCN1600, respectively, from the CCN400 value. 


Insight into these differences in cloud properties can be 
gained by examining the changes in budget terms (Fig. 7). 
Cloud condensation rates increase, but not enough to 
compensate for the decrease in cloud nucleation rates, 
such that there is less production of cloud water with 
increasing cloud-nucleating aerosol concentrations. The 
trends in these two terms oppose one another since nu- 
cleation and condensation compete for the available 
water vapor. In the more polluted cases, total condensa- 
tional growth increases because of the increased total 
surface area of the more numerous but smaller cloud 
droplets (not shown). This change effectively decreases 
the supersaturation and inhibits cloud nucleation. The 
overall increase in cloud water mixing ratio seen pre- 
viously is ultimately a result of decreased riming caused 
by the decreased size of the cloud droplets. 

There are many similarities between warm-frontal and 
orographic clouds — a sloped surface, weak to moderate 
updraft speeds and the presence of a mixed phase to 
name three — such that a comparison with results from 
orographic cloud sensitivity studies is useful. The trends 
in the three terms discussed above are consistent with the 
results of Saleeby et al. (2009) who looked at the impacts 
of CCN on an orographic snowfall case. In their study 
though, they found decreased sensitivity to riming rates 
for CCN concentrations above 500 cm~ 3 , whereas in the 
present study significant decreases in riming are seen 
even for the CCN1600 case (Fig. 7). 

Interestingly, there is very little change in cloud-to- 
rain processes across the three simulations such that 
they do not contribute to changes in mixing ratio of 
cloud water or rain. The same result was seen in one of 
the simulations of an orographic mixed-phase rain case 
examined by Muhlbauer et al. (2010). However, this 
result is in contrast with other purely warm rain studies 


that show suppressed rain formation in polluted sce- 
narios (e.g., Albrecht 1989; Xue and Feingold 2006; 
Saleeby et al. 2010). Figure 8 shows average vertical 
profiles of the rain mass produced through cloud-cloud 
collisions (autoconversion, Fig. 8a) and the rain mass 
growth through cloud-rain collisions (accretion, Fig. 8b). 
It can be seen that the autoconversion of cloud droplets 
to rain decreases as expected because of the smaller cloud 
droplet size, but that the accretion of cloud droplets by 
rain increases with increasing aerosol concentrations. 
Since formation of raindrops in this case occurs primarily 
through melting, they can gain large sizes without grow- 
ing through autoconversion of cloud water. Collision ef- 
ficiency between cloud droplets and raindrops is not 
usually very sensitive to cloud droplet size, so rain drops 
will be able to accrete more cloud water in the more 
polluted cases despite decreases in cloud droplet size 
since more cloud water mass is available because of the 
decreased autoconversion. It is unclear whether the 
nearly exact cancellation between cloud droplet auto- 
conversion and cloud-rain collisions in this case is 
fortuitous or not. Additionally, note that about twice as 
much rain mass is gained through cloud-rain collisions 
than through autoconversion. This reemphasizes the 
fact that very little rain is being produced independent 
of ice processes. 

In the more polluted simulations PSA mixing ratio 
increases while that of GH decreases, as discussed above. 
The budget terms in Fig. 7 show that vapor deposition 
increases, while riming decreases with increasing aerosol 
concentrations. Vapor deposition occurs primarily on the 
PSA species, so the increase in that process alone can 
explain the increase in PSA mass. Most of the change in 
vapor deposition is confined to the mixed-phase region of 
the cloud (Fig. 8c). This change in vapor deposition may 
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Fig. 8. Average vertical profiles of (a) rain mass growth through cloud-cloud collisions (autoconversion), (b) rain mass growth through 
cloud-rain collisions (accretion), and (c) ice mass growth and decay through vapor deposition and sublimation. 


be due to an alteration in the competition for vapor be- 
tween liquid and ice particles in this region because of 
changes in the size and number of cloud droplets. In ad- 
dition, riming is one pathway for the conversion of PSA 
to GH; therefore, the decrease in riming further explains 
the increase of PSA mass as well as the decrease in GH 
mass with increasing aerosol concentrations above the 
melting level. At and below the melting level, remaining 
PSA mass is quickly converted to GH and then to rain 
(Walko et al. 1995). The changes in vapor deposition and 
riming almost exactly cancel one another, so below the 
melting level there are no significant differences in the 
surviving ice mass or in the rain mass produced through 
melting. 

The two largest contributors to rain mass, melting 
and collision-coalescence, both show little sensitivity to 
cloud-nucleating aerosol concentration, as discussed 
above. This explains why there is also little change in 
rain mixing ratio below the melting level across the three 
simulations. As would be expected in association with 
this result, there are also only small changes in the av- 
erage precipitation rates (Table 1). 

b. Local precipitation changes 

Although on average the precipitation rates are not 
very sensitive to aerosol concentration, there are larger 
changes seen locally. Fig. 9a shows the time averaged 
precipitation rates as a function of distance from the 
front. The largest changes are seen between 50-200 km 
from the front where there is on average a 9.5% differ- 
ence between CCN400 and CCN1600; however, the 
trends are often nonmonotonic and it is difficult to at- 
tribute any specific processes to the changes. Therefore 
we would not expect the location of these changes to be 
a general result of this study. 


Figure 9b again shows time-averaged precipitation 
rates but excludes those rates greater than 5 mm h 1 
that may be associated with more convective-type pro- 
cesses. The local maxima in CCN400 and CCN1600 in 
Fig. 9a around 50 km from the front are effectively re- 
moved and much clearer trends in precipitation are re- 
vealed. There is a delay in the onset of precipitation with 
increasing cloud-nucleating aerosol concentration such 
that precipitation is initially shifted about 15-20 km 
downwind of the front. This is similar to the “spillover 
effect” seen in many orographic snow cases in which 
precipitation is shifted to the leeward side of a mountain 
in more polluted conditions (e.g., Givati and Rosenfeld 
2004; Jirak and Cotton 2006; Lynn et al. 2007; Saleeby 
et al. 2009). However, farther from the front rain rates in 
the two more polluted simulations have higher in- 
tensities at their peaks. Just after these peaks, rain rates 
are similar across all three cases. Because of the changes 
in peak rain rates, the difference in average rainfall is 
still only 2.0% between CCN400 and CCN800 and 1.4% 
between CCN400 and CCN1600. These are not signifi- 
cant changes in precipitation. This is in contrast to some 
orographic sensitivity studies which generally find de- 
creases in the overall precipitation because of losses 
through sublimation of snow and ice crystals on the lee 
side of the mountain (Lynn et al. 2007). However, in this 
case, since there is little precipitation where sublimation 


Table 1. Summary statistics of precipitation rate. 



Precipitation rate 

Change from 

Simulation 

(mm h _1 ) 

CCN400 (%) 

CCN400 

1.041 

— 

CCN800 

1.056 

1.4 

CCN1600 

1.06 

2.1 
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Fig. 9. (a) The average precipitation rate as a function of distance from the front, (b) As in (a), but excluding 

precipitation rates greater than 5 mm h 1 . 


of ice and evaporation of cloud water is substantial, 
precipitation cannot be reduced significantly through 
these processes. 

Examination of the vertically integrated processes 
contributing to rain formation in areas where the rain rate 
is 5 mm h 1 or less (Fig. 10) shows that both cloud water 
conversion to rain and melting are suppressed near the 
front. The suppression from melting is not seen between 
the CCN800 and CCN1600 cases, but suppression is still 
clearly evident in the cloud-to-rain processes. It is not 
immediately obvious why there is so little difference in 
melting between CCN800 and CCN1600, especially since 
riming does not show the same insensitivity. However, 
the GH particles are larger in CCN1600 and their faster 
sedimentation rates may be off setting the decrease 
in mass because of decreases in riming rates. The in- 
creased peak in precipitation in CCN800 and CCN1600 
farther downwind appears to be because of enhance- 
ments in both melting and cloud-to-rain processes at 
these locations. 

c. Latent heating and warm-frontal structure 

The total latent heating (Fig. 11) shows only small 
differences among the experiments, despite the fact that 
significant changes were seen in the rates of many cloud 
processes. There are increases up to 1 km, and again 
between about 3 and 3.5 km. The increase at this latter 
level is about 45% from CCN400 to CCN1600, which 
is quite significant, although very localized. Decreases 
are seen around the melting level as well as near 4 km. 
When vertically integrated, the difference between 
CCN400 and CCN1600 is only 1.5%. This implies that 
though changes in warm-frontal latent heating from 
specific cloud processes vary significantly, they largely 
cancel one another when summed. 


Finally, vertical profiles of vertical velocities are shown 
in Fig. 12. There is on average a small (3%-5%) increase 
in vertical velocity throughout the column with increasing 
aerosol concentrations. Increases in updraft speed (Fig. 
12b) are only seen in CCN1600 and only between about 
4 and 6 km, possibly because of the localized increase in 



Melting 



Distance from Front (km) 


Fig. 10. Vertically integrated and averaged values of selected 
process rates as a function of distance from the front. 
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Fig. 11. Average vertical profiles of the total latent heating rate. 


latent heating seen from 3 to 3.5 km. Also note that the 
level of maximum latent heating is below that of maxi- 
mum vertical velocity. This may be due to lower vapor 
availability where updrafts are stronger aloft. Downdrafts 
(Fig. 12c) show more variability and clear, monotonic 
trends are not evident. Overall, the changes to average 
vertical velocity are small and do not seem to have had 
significant impacts on the warm front as a whole. 

6. Conclusions 

RAMS was successfully used to simulate the impacts 
of increased aerosol concentrations on a springtime 
extratropical cyclone that developed and crossed the 
United States from 9 to 11 April 2009. Although 


we simulated an actual case, the results of the aerosol 
sensitivity tests can be treated in an idealized framework. 
First, the budgets of cloud ice, cloud water, and rain were 
examined over the warm front. Vapor deposition and 
riming both contribute significantly to the total ice mass 
in this case since updraft speeds were high enough to 
sustain a deep mixed-phase cloud region. Riming is the 
most significant sink of cloud water, especially in the 
mixed-phase region. The melting of ice species produces 
75% of the rain with the other 25% coming from con- 
version of cloud droplets to rain. The exact contribution 
from each of these processes varies as a function of dis- 
tance from the front, as seen in Fig. 10, with conversion of 
cloud droplets being fairly constant from 25 to 175 km 
ahead of the front and melting peaking sharply at 200 km 
ahead of the front. These results, while specific to the case 
analyzed here, do generally agree with the studies of 
Rutledge and Hobbs (1983) and Gedzelman and Arnold 
(1993) who both also examined cases in which rain was 
the predominant hydrometeor reaching the surface. This 
suggests that the results presented here are robust for 
warm fronts in which the seeder-feeder process is dom- 
inant. The budgets would be expected to differ for colder 
systems that do not produce rain, or for systems in which 
convection is more prominent. 

Three sensitivity experiments were conducted in 
which the profiles of cloud-nucleating aerosol concen- 
trations decreased linearly to 100 cm~ 3 at 4 km AGL 
starting with surface concentrations of 400, 800, and 
1600 cm~ 3 . There is a slight shift seen in the distribution 
of light to moderate precipitation rates where less pre- 
cipitation falls with increasing aerosol concentration at 
the surface position of the front, but more falls between 
100 and 200 km from the front. In total there is a mere 
2% increase in the temporally and spatially averaged 
warm-frontal precipitation from CCN400 to CCN1600. 

Significant changes in association with variations in 
aerosol concentrations are seen in many processes that 
contribute to the production and depletion of cloud 
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Fig. 12. Average vertical profiles of (a) total vertical velocity, (b) updrafts, and (c) downdrafts. 
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Fig. 13. Summary schematic showing the changes to the microphysical properties of the 
mixed-phase portion of the warm-frontal cloud and the total precipitation. Vapor deposition 
onto pristine ice, snow, and aggregates (PSA) increases (as indicated by the change in arrow 
thickness) while riming of the mixed-phase species (GH) decreases. These two trends in the 
growth of ice mass cancel one another resulting in little change in total ice mass and little 
change in rain production through melting. Differences in size and number of hydrometeors 
between the clean and polluted scenarios indicate qualitative changes. 


water and ice, but not of rain. A simple schematic is 
shown in Fig. 13 that summarizes these results. Because 
of higher cloud-nucleating aerosol concentrations, cloud 
droplets are smaller and more numerous. In the mixed- 
phase region of the cloud, these changes in the cloud 
droplet distribution appear to cause enhanced deposi- 
tion of vapor onto PSA crystals but decrease riming 
by GH species. The responses of vapor deposition and 
riming tend to cancel one another such that similar 
amounts of total ice are produced, causing little change 
in rain production from melting, which is the most 
important source of rain in warm-frontal clouds. This 
suggests that the mixed-phase part of the cloud has 
compensating responses that limit the impact of increased 
cloud-nucleating aerosols on surface precipitation. In 
addition, the compensating responses in cloud pro- 
cesses prevent large changes in total latent heating. As 
a result, cloud-nucleating aerosols do not appear to be 
able to play a strong role in determining the total pre- 
cipitation, average latent heating, or average vertical 
velocities of warm fronts. 

Only one case of a spring-season, mostly stratiform 
warm front was examined in the present study, and it is 
unclear whether these results can be generalized to all 
warm fronts. The response of more convective fronts or 
colder systems with less liquid water may be different, as 
may the response to ice nuclei or giant CCN for storms 
over the ocean. Investigation into these questions is left 
for future work. However, an increasing number of sen- 
sitivity studies (Stevens and Feingold 2009; Muhlbauer 
et al. 2010; van den Heever et al. 2011; Seifert et al. 2012), 
including the one here, are finding that the mixed phase 


of cloud systems provide built-in buffers that resist 
wide-scale changes in aerosol concentrations. This trend 
lends support to the robustness of the conclusions pre- 
sented here. 
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